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Abstract. The latest generation of volatility derivatives goes beyond vari- 
ance and volatility swaps and probes our ability to price realized variance and 
sojourn times along bridges for the underlying stock price process. In this pa- 
per, we give an operator algebraic treatment of this problem based on Dyson 
expansions and moment methods and discuss applications to exotic volatility 
derivatives. The methods arc quite flexible and allow for a specification of 
the underlying process which is scmi-paramctric or even non-parametric, in- 
cluding statc-dcpcndcnt local volatility, jumps, stochastic volatility and regime 
switching. We find that volatility derivatives are particularly well suited to be 
treated with moment methods, whereby one extrapolates the distribution of 
the relevant path functionals on the basis of a few moments. We consider a 
number of exotics such as variance knockouts, conditional corridor variance 
swaps, gamma swaps and variance swaptions and give valuation formulas in 
detail. 



1. Introduction 

Volatility derivatives are designed to facilitate the trading of volatility, thus 
allowing one to directly take a range of tailored views. A basic contract is the 
variance swap which upon expiry pays the difference between a standard historical 
estimate of daily return variance and a fixed rate determined at inception, see 
(Dupire 1992), (Carr and Madan 1998) and (Derman et al. 1999). A variant on 
this is the corridor variance swap which differs from the standard variance swap 
only in that the underlyings price must be inside a specified corridor in order for 
its squared return to be included in the floating part of the variance swap payout 
(Carr and Lewis 2004). A further generalization is the conditional variance swap 
which pays the realized variance of an asset again within some corridor, whereby the 
average is taken only over the period when the spot is in the range. The advantage 
of conditional corridor variance swaps is that they allow one to take a view on 
volatility that is contingent upon the price level, gaining exposure to volatility only 
where required. Although conditional variance swaps appear to be traded, see (J. P. 
Morgan Securities 2006), this is the first article in the open literature proposing a 
pricing methodology. 

The numerical method we present extends without difficulties to other exotic 
volatility contracts such as variance swaptions, see also (Carr and Lee 2007). We 
also take the liberty of inventing exotic volatility derivatives such as variance knock- 
out options. These do not seem to be much traded although perhaps they should be. 
A variance knockout can be regarded as a variation on the theme of barrier knock- 
out options whereby the knock-out condition is not triggered by the underlying 
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crossing a certain level. Instead, a variance knockout vanishes in case realized 
variance prior to maturity exceeds a certain pre-assigned threshold. The benefit of 
a variance knockout over a barrier knockout is that its hedge ratios are smoother. 

In this paper, we use operator methods for pricing. Our presentation is self- 
contained for the specific purpose at hand but see the review article (Albanese 2006) 
for a more extended discussion of operator methods. One point we should stress is 
that the mathematical and numerical methods we present would work as efficiently 
no matter what underlying model for the stock price process is chosen. We select one 
for the purpose of discussing a concrete case and generating sample graphs, but we 
are confident the reader can do better if she intends to refine it. Our mathematical 
and numerical methods are model agnostic as they do not rely on closed form 
solvability and their performance is not linked to the model definition. Any model 
for the stock price dynamics would work just as well, the only limitation being that 
market models cannot be accommodated within the formalism we propose. 

We believe that it is important to embed econometric estimates in the volatility 
process for the underlying stock price and thus make it to reproduce realistically the 
features of the historical process. Operator methods reviewed in (Albanese 2006) 
are useful in this respect as they allow one to construct models semi-parametrically 
or even non-parametrically while resting assured that numerical efficiency is not 
affected by model choice. Our working example is a 3-f actor equity model which 
encompasses stochastic outlook dynamics, stochastic volatility, local volatility and 
jumps of both finite and infinite activity. This process can formally be expressed 
as follows: 

(1) dSt = fit{St)dt + c7t{St)dWt + jumps. 

Here both the drift and the volatility terms can be freely specified, thus making 
available additional degrees of freedom and allowing for a time-homogeneous cali- 
bration. 

Armed with the calibrated model we proceed to a moment method which allows 
us to obtain conditional moments of integrals of stochastic processes. Having the 
first few moments as a starting point, the distribution is extrapolated by matching 
moments with elementary probability distributions functions known in analytically 
closed form. On this basis we then price volatility derivatives and find their hedge 
ratios. 

The article is organized as follows: The next section describes the underlying 
model, and the following one the moment method. Sections Q, ([s]) and ^ de- 
scribe the variance swaps, corridor variance swaps and conditional variance swaps 
respectively, and their pricing within this framework. In section ([t]) we look at 
some more variance related contracts, particularly the gamma swap and variance 
knock-out options. A final section concludes. 

2. Underlying Model 

Our working example, as a base model, is specified semi-parametrically and 
amounts to a 3-factor equity model with stochastic volatility and stochastic outlook 
dynamics. We find that postulating a slowly varying stochastic outlook dynamics 
alongside a faster mean-reverting volatility dynamics is essential to mimic the his- 
torical features of the process and calibrate to both short-dated and long dated 
forward volatiHty. The presence of jumps helps to explain the shorter maturities, 
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stochastic volatility for the medium to long term and the outlook dynamics for the 
long-term features of the model. 

Near time-homogeneity as a key property and requirement for our model; apart 
from an exogenously specified term structure of interest rates, we achieve high 
precision calibration to all of our objectives within the bounds of a nearly time- 
homogenous model. We stress that time-homogeneity is not a technical constraint 
from the numerical point of view, but rather an economic constraint to discipline 
the calibration of otherwise very flexible model specifications. Since parameter 
sensitivities and hedge ratios one obtains numerically are remarkably smooth even 
if one uses single precision floating point arithmetics on GPUs, the flt to calibration 
targets can be made arbitrarily accurate by reflning the model with modest time- 
inhomogeneities . 

The model is specifled in terms of a continuous-time Markov chain on a lattice 
with flnitely many states. We select the following notations for our discretization: 
the base lattice A is made up by triples y = (x, a, 6) where x = . . . 69, a takes 
values stable and negative is an indicator of outlook regime and h takes values 
low, medium and high is an indicator of volatility regime. The stock price is given by 
a function S{x). See (Albanese 2006) for the theory behind these discretizations, in 
particular for rigorous convergence proofs and sharp bounds on convergence rates. 

We specify the Markov generator to precisely express the model sketched in 
equation ([T]). This matrix then needs to be exponentiated numerically to obtain 
the required transition probability kernels. We are therefore describing a three- 
dimensional problem: one for the underlying stock price, one for the stochastic 
volatility regime and one for the stochastic drift regime. Since we have multiple 
outlook regimes characterized by different drifts, we match the jump amplitudes 
non-parametrically to equal the risk- neutral drift. We approximate interest rates 
r{t) as deterministic and piecewise constant in such a way to match the spot dis- 
count function. With the interest rates taken to be piecewise constant over time 
intervals of the order of a quarter, we specify a Markov generator Ci on each such 
interval [ti,ti^i) as having the following form: 

£(x, a, 6; x' , a', b';t) = {Vontix; x'\a, h)5aa' - Joyit{x, a; x' , a'\h)) Sbb' 

(2) -I- Jvg{x; x'\b)Saa'Sbb' + Vsv(6; b'\x, a)6xx'Saa' 

The tri-diagonal matrix I?out(a:; x'\a, b) is chosen so that we match the risk neutral 
drift and satisfy the no-arbitrage condition in equation ^ below. The operator 
'Jouti^T o,; x' , a'\b) models jumps downwards of exponential type which in our ex- 
ample, are chosen to be of size 12% if the outlook is negative and small symmetric 
jumps of size 2% if the outlook is stable. Both jumps are flnite activity and occur 
with yearly frequency. 

The operator J'vg{x;x'\b) describes a variance-gamma (inflnite activity) process 
of log-normal volatility dependent on the state variable b and of the form S^^^''~^ 
where P{S) is the function in Fig. [l] The variance rate is chosen to be 4% and is of 
the order of the rate estimated econometrically, not sufficient to sustain an implied 
volatility skew over medium or long time horizons. The first two terms admit the 
state a; = as an absorbing point. 

The operator Vsv{b;b'\x,a) describes a mean reverting process for stochastic 
volatility given by volatility regimes. It is allowed to switch between three possible 
states: low, medium and high. Evidence that the volatility process experiences 
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patters which can be effectively described by means of regimes has been previously 
discussed by Rebonato and Joshi in (Rebonato and Joshi 2001). The underlying 
stock price is allowed to transition from one regime to another only between adjacent 
regimes. 

The Markov generator satisfies the following three conditions corresponding to 
probability positivity and conservation and to arbitrage freedom, respectively: 

(3) C{x,a,b;x',a',b';t)>0 if {x,a,b) ^ {x' ,a' ,b') 

(4) ^ C{x,a,b;x',a',b';t) =0 

(5) C{x,a,b;x',a',b';t){S{x')-S{x))=r{t)S{x) 

a:' ,a' 

The problem of computing the propagator U{t\,t2) satisfying the backward equa- 
tion 

(6) ^U{h,t2) + C{h)U{h,t2)=i) 

at I 

with final time condition U{ti,t2) = I, the identity operator, is solved by exponen- 
tiating our generator jC{t) 

(7) U{ti,t2)=exp(^JJ £{s)ds^ 

The matrix exponentiation can be done in a number of ways (Moler and Loan 2003) . 
The method we advocate is the linear fast exponentiation algorithm which we find 
is fast, stable with respect to floating point errors and can be implemented in single 
precision on highly efficient GPU platforms. 

Fast exponentiation proceeds as follows: Assume that the generators £{t) are 
piecewise constant as a function of time. Suppose jC{t) = Ci in the time interval 
\ti,ti -\- (Af)i]. Assume 5t is chosen to be small enough (but not smaller) that the 
following two conditions hold: 

(8) imiiy^^{l + 5tCi{y,y))>l/2 




= n e N 
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This condition leads to intervals St of the order of one hour of calendar time. To 
compute e^^^^^^^{x,y), we first define the elementary propagator 

(10) ust{x,y) = d^y + dt£i{y,y) 
and then evaluate in sequence 

(11) U2St = Ust ■ Ust , UiSt = U25t ■ U2St , • ■ • , = U2"~i5t ' U2"-^St 

The matrix multiplication is accomplished numerically by invoking the level-3 BLAS 
routine sgemm. 



3. The Moment Method 

One possible strategy to evaluate moments of path dependent functions is to use 
Dyson's formula, see (Albanese 2006). This formula is sometimes also attributed 
to Kac (Darling and Kac Mar., 1957) and provides an expression for moments of 
integrals of stochastic processes. Moment methods involve reconstructing probabil- 
ity distribution functions for such integrals by extrapolating out of the knowledge 
of the first few moments. This involves selecting a probability distribution function 
available in closed form and matching moments. 

Consider a time interval [T,t\ and a Markov generator h(yi,y2',t), and consider 
the integral It given by 

(12) It = I{yut)= f ^iys,s)ds 

Jt 

Let's introduce the following one parameter family of deformed Markov operators 
parameterized by e e M 

(13) L^{yi,y2;t) = L{yi,y2;t) + eV{yi,y2;t) 
where 



(14) V{yi,y2;t) = 4>{yi\t)5y^^y2 + L{yi,y2;t)x{yi,y2;t){ij{y2;t) - ip{yi\t)) 



The functions 4){y; t), ipiy] t) ^iid xiV) ^) ^re all continuous functions. We add the ^ 
and X functions here as we can use them to model certain path dependents, however 
in the applications we describe below we only need the first term (j)(y; t) in order to 
model the integral in equation (12 1. 

Theorem 1. For n = 0, 1, 2, . . . , we have that 



(15) 



-Pexp (^j^ L,{s)ds^{yi,y2) = ET[l^S{yt - y2)\yT = ■ 



Proof. Consider Neper's formula for the propagator 

N 



(16) 



Pexp(/^M.)d.)=^lim Pn(l 



t-T 
N 



(L(t,) + ecp{ti)) 



N 
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where ti = T + ^ . By collecting similar powers in e, one finds Dyson's formula 
(17) 

Pexp f j h^{s)ds] = exp((f - T)L+ 

(18) 



(19) E^" / rf,si... / ds„fe(^^-^)^F(si)e(«=-«^)^....y(s„)e(*- 
The time-ordered integrals above are proportional to conditional moments, i.e. 

(20) Pexp L,(s)dsj (yi, 2/2) = X] ^^^^ [/("^(j/t " 2/2) |yT = yi] • 

Here, the factorials originate from the time ordering. □ 

A technique which is numerically stable in many situations of practical rele- 
vance is to numerically differentiate with respect to e the deformed propagators 

Pexp ^ Le(s)(is^ (2/1, 1/2) and evaluate at e = 0. This technique can be used 

to obtain the first moments of It on any given bridge for the underlying Markov 

process. In most applications, wc find that two moments suffice to extrapolate the 
probability distribution function to sufficient accuracy. To do so, it is convenient 
to choose from among the probability distribution functions which are analytically 
tractable. For instance, starting from the first two moments only, one can use cither 
the log-normal or the chi-square distribution. In the above theorem, one only needs 
to compute the formula for n equal to the number of moments we want to use. 

One may notice that Dyson formula in our Finance application is used differently 
from how it is used in Physics. In Physics one usually has a perturbative method 
to evaluate moments and the purpose is to reconstruct the kernel. In Finance, fast 
exponentiation gives a good way to evaluate the kernel and Dyson formula is used 
to evaluate moments by numerical differentiation. 



4. Contracts on Realized Variance 

Having introduced a stock price model and a moment formula for path function- 
al, we next consider volatility derivatives. 

A variance swap is a financ;ial coiitrac;t that upon expiry pays the difference 
between a capped measure of the historically realized variance of daily stock returns 
and a fixed swap rate determined at inception. As usual, the swap rate is initially 
struck at equilibrium so that the variance swap has zero cost to enter (Carr and 
Madan 1998). 

Variance swaps of maturity t and time of issuance T have a payoff given by 



(21) min y* v{ys, s)ds, f ■ SR''^ - SH 
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where SR is the swap rate and / is a factor, a typical value being / = 6.2. Here, 
v{yi,t) is the instantaneous variance defined as follows: 

(22) v{y,,t)=J^C{y^,y2;t)log'(^^^ 

The cap is customary on single stocks as, in case of default, realized log-normal vari- 
ance may actually diverge. Instead, variance swaps on indices are usually uncapped 
and the equilibrium swap rate is simply given by the first moment 



(23) —E, 



v{ys, s)dsd{yt - 2/2) 



yr = yi 



SR' 



The expectation is given by the first derivative in equation ( 15 1 without any need for 
moment matching. Summing over all final points ?/2 removes the bridge condition 
on the final points. 

The payoff of a volatility swap is 



^-3^ J v{ys,s)ds,SRj - SR 



To price these contracts one needs to evaluate the distribution of realized variance 
on a bridge, i.e. of the functional 

(24) RV{y2) = 5{yt-y2)j^ j^v{ys,s)ds 

In this case, moment matching is required even in case the cap is lifted. 
4.1. Matching the moments to a chi-square distribution. Let 

mi = Et [it 5{yt-y2)\yT = yi] and m2 = Et [it S(yt - y2)\yT = yi] 

be our first and second pre-assigned moments. The standard chi-square distribution 
is given by 

■'^ ' 2r(f) V2/ 

where a is the number of degrees of freedom. The first and second moments of this 
distribution are 

E[x\^a, E[x']^a{a + 2) 

To match two pre-assigned moments mi, m2i one can pass to the new variable 

mi 

C = — X 
a 



and chose 



2m\ 



m2 — m\ 



The cumulative distribution function is 

where T{z) and 7(2, a) are the gamma and incomplete gamma functions respec- 
tively, i.e: 

/■oo pa 

T{z)^ s^'-^e-'ds, ^{z,a)= s'^'^e-'ds. 
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We find that 



Et [min(W, iJKnax) 5{yt - 2/2)] = ^ - F{K- a)) + ai^(if; a + 2)] 



and 



Et [min(x/i?F, /RIW) ^(yt - ^2) 
where 

a = a(yi,2/2) 



if (1 - a)) 



2TOi(2;i,y2)' 



1(^1,2/2) 



. N / \2' K=K{yi,y2)= Wmax 

m2{yi,y2j - mi{yi,y2)^ wi (2/1,^2) 

These formulas allow one to price both variance swaps and volatility swaps. Since 
the dependency on the swap rate in both cases is non-linear, an exact calculation 
requires using a root finder. This may be avoided as an excellent approximation is 
often obtained by estimating the variance swap rate as follows: 



(25) 



VARSRt = ET[ram{RV, i?Knax) S{yt - 2/2)] 



and the volatility swap rate as follows: 

(26) VOLSRt = J2^t UmiVRV, y^RV^) 5{yt - 2/2) 



where 
(27) 



2/2 



-RVmax = fYU{yT,T; yz,t)'m2{yT,yz) 
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Figure 2. Volatility and variance swaps (Chi-square) 
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4.2. Matching the moments to a log-normal distribution. Let x be log- 
normally distributed with probability distribution function 



1 



-(ln2;-p)^/2(T^ 



Xa^/2lT 

then the first two moments are given by 

Given these expected values, we can obtain relationships for /i and a as 



(28) 



^l = log 



E\x\ 



and cr — log 



E[xY 



and take the pre-assigned moments to be equal to these expected values E[x\ = wi 
and Elx"^] = 1712- We obtain 

ET[mm{RV, i?K„ax) S{yt - j/a)] 

expressed in terms of the well known standard normal cumulative distribution func- 
tion as 

(29) e''+^V2_^ j^ log(flK.ax)^-M-^V2 j ^ ^^^^^^^^ |^ log(W^..)-M j 
and the second expectation is 



Ej 



M^RV, v/i?K„ax) ^(yt - 2/2) 



(30) ei(^^+-^)AA (^ l°g(^^"-)^- - -^V2 ^| ^ yj^^ j^ log(i^Knax) 
where /i and cr are specified above in terms of mi and m2. 
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Figure 3. Volatility and variance swaps (log-normal) 
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4.3. Matching higher moments with the Pearson distributions. The log- 
normal and the chi-square distributions allow us to match the first two moments, 
here we look at three. Consider the Pearson Type III distribution which has prob- 
ability distribution function 

defined on [a, -l-oo). The special case of this, when a = 0, h = 2 and p is half of an 
integer, gives the Chi-Squared distribution. In general, the moments are given by 

(31) E[x\ ^a + bp 

E[x^] = {a + bpf + b^p 

E[x^] = {a + bpf + Sb'^pia + bp) + 2b^p 

and matching these with our pre-assigned moments mi, m2 and ma and computing 
the values of a, 6 and p in terms of these moments we find 

2(m2 — ml)^ 



(32) a = nil 



ma + 2m\ — 3mim2 



^ ma -|- 2m^ — 3mim2 
2(m2 — m\) 
4(m2 — m?)^ 
(ma -t- 2mJ — 3mim2)^ 
Then the expectation in concern, expressed in terms of the chi-squared cumulative 



distribution function F{x;a) (section 4.1) is 

RV, 



Et [Tnin{RV, i?K.ax) 5{yt - y2)] = {a + bp- i?Knax)i^chi 2p, 2 



max 



(33) ^bi^-^ j ^^^^ + (2i?Knax -a- bp) 

5. Corridor Variance Swaps 

Corridor variance swaps differ from standard variance swaps in that the under- 
lyings price must be inside a specified corridor in order for its squared return to be 
included in the fioating part of the variance swap payout (Carr and Lewis 2004) . 

Corridor variance swap are sometimes engineered so that they afford full protec- 
tion against situations where the variance spikes accompanied by a decline in stock 
value, but offer only limited protection in case the volatility surge is associated to 
a stock price rally. Following (Carr and Lewis 2004), one could achieve a robust 
static replication by means of a strip of options as in the case of the variance swap, 
only that here the hedge is limited to the pre-specified corridor. 
The payoff of corridor variance swaps are given by 

1 /■* 

(34) — v{y,,s) liL < Sivs) < H)ds - SR^ 
so we are interested in the integral 

1 /■* 

It = -j-—^ / (t>{ys,s)ds 
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where 
(35) 



y2 



is the "instantaneous corridor variance", and solve the problem of computing its 
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Figure 4. Corridor Variance Swaps 

moments by specifying in this case the operator 

^e{yi,y2;t) = C{yi,y2;t) + e(t){yi,t)6y^y^. 

From here the methodology used in the case of ordinary variance swaps carries 
through and one is reduced to computing and matching moments. 

6. Conditional Variance Swaps 

Conditional corridor variance swaps are an interesting variant as they can more 
easily be interpreted in intuitive terms. In this case, the payoff is given by the 
accrued variance divided by the time elapsed in the range, not the total number of 
days (J. P. Morgan Securities 2006). Namely, the payoff is given by 



(36) 



Jlviy,,s)liL<S{ys) <H)ds 
J^l{L<S{ys) <H)ds 



SR' 



The principal use of conditional derivatives is when an investor expects a particular 
asset (e.g. market level, volatility) to outperform contingent on a level of another 
asset being breached. In the case of a conditional variance swap, the investor may be 
anticipating that volatility will increase/decrease if the market were to rally/sell-off 
or stay within a range (J. P. Morgan Securities 2006). While an investor can select 
the corridor freely, two examples are the conditional up- variance and down- variance 
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swaps. Up-variance accrues realised variance only when the underlying is above a 
pre-specified level (i.e. no upper barrier), while down- variance is accrued only when 
the underlying is below the specified barrier (i.e. no lower barrier). 

To apply the moment method to this case, the first thing to note is that essen- 
tially we are modelling the two integrals appearing in the payoff at the same time. 
We are going to need a bi-variate distribution to do this. Let's write: 

^ = £ vivs, s) 1(L < S{ys) < H)ds 



r(2) 



= j 1{L< S{ys) < H)ds 



and in order to compute the expectation 

(37) E 

we'll need the following expectations: 
E 



r(2) 







V(2)' 


, E 




, E 




, and E 


7(1) r(2)' 
-'t -'t 



To tackle this problem consider the operator 
(38) 
where 



(39) 

and 

(40) 



(2/1,1/2) = C{yi,y2) + ei<P{yi)Sy,y^ + £2 V' (2/1) ^3/13/2 
cf>{yut) = ^/:(2/i,y2;i)log^ (lllj) 1{L < S(yi) < H) 



V'(yi) = 1{L < S{yi) < H) 

then the problem of computing these moments is solved using the above procedure. 
In particular, we have 

d 



aei 



e«*-^)^--)(2/i,2/2)=i? 



ei=0 



r(i) 



yr = yi,yt = 2/2 



and 



del 



ei=0 



similarly (but with respect to €2) for 



E 



(2) 



yr = yi,yi = y2 



and E 



(4'r 



vt = yi,yt = 2/2 



yr = yi,yt = y2 



The joint expectation will involve the mixed derivative 



(41) 



E 



r(l) r(2) 
-'t -'t 



d€ld€2 



e«*-^)^--)(2/i,2/2) 

«i,e2=0 

and once computed, we make use of these expectations to match a bivariate dis- 
tribution. For simplicity of notation we leave out the conditional part of these 
expectations, noting that all the moments we obtain axe conditional to the initial 
and final points. 
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6.1. Bi-Variate log- normal distribution: Let 

Yi \ / log(Xi) 

Y2 J { log(X2) 

then Xi and X2 have joint probabihty distribution function (Johnson and Kotz 
1942): 

1 



(42) 



exp ■ 



2(1 



logxi - ^1 

0-1 



2TTaia2\/l - p'^xiX2 



-2p 



{\0gXi- lJLi){\ogX2- 1^2) , flogX2~fJ-2 



(TiCr2 



0-2 



the bivariate log-normal distribution, where p is the correlation between Xi and 
X2. Both and X2 are log-normally distributed with 



E[X,] = 



and E[Xf] = 



2] _ 2Aii+2cr^ 



i = 1,2. 



Matching these with the pre-assigned moments, and solving for fii and Cj we find 



(similarly to section 4.2) 
2 / / 



= log £; 



(0 



and erf = log 



Moreover, the mixed term is 

E [X1X2] E [e^'+^^] = e^i+^^+3('"i+'"'+2'"^i'"^) = E[Xi]E[X2] 6^"'" 
which gives p (in terms of the pre-assigned moments): 



(43) 



P 



aiCT2 



log 



E 


"r(l)r(2)' 




E 




E 


r(2) 



Now we are in a position to compute the expectation E 
bivariate log-normal distribution, we have 



,(1) / ,(2) 



Given the 
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Xn 



= i?[e^i-^^] ^E 



where K 



2 — ■~Y2 will also be normally distributed ~ N{—p2, 0-2)1 
Xi 



E 



X2 



and the expectation ( 37 1 is given by 

E (l 
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Plugging this in the payoff of the conditional variance swap ( 36 1 gives us the value 
of the conditional variance swap. 

We could go a step further and compute the analogy of the capped variance swap 



in equation (21 1 where we are now interested in the payoff E 
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Figure 5. Conditional Variance Swaps: A standard variance 
swap, two conditional variance swaps of different corridors, and 
a f 05% (of spot) conditional up-variance swap. 



Writing this expectation as a double integral over the bi-variate probability distri- 



bution function (42 I we find 
(45) 



E 



X. 



Ma 



X. 



Ma 



This double integral will need to be evaluated numerically; one could apply a Gauss- 
ian quadrature. 



7. Gamma Swaps and Variance Knockouts 

In this section, we apply the moment method to other volatility contracts, namely 
options on realized variance, variance knock-out options and gamma swaps. 

7.1. Options on Realized Variance. Let's consider the example of a call option 
on realized variance, we are interested in the payout E — A')+]. Computing 
the conditional moments for the realized variance as described above, we can sum 
over all final points to obtain the unconditional moments. Matching these with 
a distribution we know, as done above, we can look at the payoff 

(46) Et {{RY - K)-^\ 



which under the log-normal case is given by the Black-Scholes formula, with the 
correct inputs from equation ( 28 1 . 
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7.2. Variance Knock-out Options. The payoff of a variance knockout call-option 
with volatility barrier at H is 

(47) {D{T)St ~ if )+l f / s)ds < H' 



t-T 

Variance conditional knock-out options are priced as follows: 

(48) VKOt = ^iVT, T; 2/2, *) Pr (^(^2) < H^) {D{T')S{y2, t) - K)^ 

where Pr {EViy'i) < H^), the probability that RV{y2) < is obviously the cu- 
mulative distribution function (for the distribution which we chose to model RV 
with). 

7.3. Gamma swaps. The gamma swap (also known as weighted variance swap) 
is similar to the variance swap only we have an additional spot factor. They reduce 
the weighting of the downside scenarios but introduce a market level exposure. The 
gamma variance swap pays at maturity the realized weighted variance of an asset 

252 ^— \ / Si \ Si 
log 



m \Si-iJ So 

Wc note the following features of the Gamma swap 

(i) As the spot approaches zero, we do not have a blow up as in the case of the 
variance swap: limj^^o log(a;) = 0, 

(ii) Similar to the variance swap case, one can split the gamma swap payoff into a 
European payoff and a daily hedging strategy. However in the case of the Gamma 
swap, the weights in the replicating portfolio are of order 1/K compared to those 
in the variance swap of order 

To price the gamma swap within our framework, we look at the instantaneous 
level: we no longer have the instantaneous variance, but the "instantaneous gamma" 
given by 

=i:*..«.).o.^(|g)|g 

Gamma swaps of maturity t and time of issuance T have a payoff given by 



(49) — = / g{ys,s)ds-SR 



8. Summary and extensions 

Using operator theory wc presented a 3-factor equity model which is semi- 
parametric or even non-parametric and allows for a time-homogeneous calibration. 
A moment method based on Dyson's formula is then given and used to compute 
the moments of integrals of the underlying Markov process. We then apply this 
method to a various volatihty derivatives which we find to be well suited to moment 
methods. We consider a number of exotics such as variance knockouts, conditional 
corridor variance swaps, variance swaptions and gamma swaps. Straightforward ex- 
tensions would include for instance conditional corridor variance knockout options, 
forward starting variance swaptions, etc.. 
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Figure 6. Gamma swap 
References 

Albanese, C. (2006). Operator Methods, Abelian Processes and Dynamic Conditioning, preprint, 

available at www.level3finance.com. 
Carr, P. and D. Madan (1998). Towards a theory of volatility trading. Volatility, Risk Publications, 

ed. Jarrow, R. 

Carr, P. and K. Lewis (2004). Corridor variance swaps. Risk February, 67—72. 

Carr, P. and R. Lee (2007). Realized Volatility and Variance: Options via Swaps. Risk. 

Darling, D. A. and M. Kac (Mar., 1957). On occupation times for markoff processes. Transactions 

of the American Mathematical Society 84, 444-458. 
Derman, E., K. Demeterfi, M. Kamal and J. Zou (1999). More Than You Ever Wanted to Know 

about Volatility Swaps. Journal of Derivatives. 
Dupire, B. (1992). Arbitrage Pricing with Stochastic Volatility, preprint, Socit Gnrale. 
Johnson, N.L. and S. Kotz (1942). Continuous multivariate distributions. John Wiley and Sons. 
J. P. Morgan Securities (2006). Conditional variance swaps, product note. Technical report. 
Moler, C. and C.V. Loan (2003). Nineteen dubious ways to compute the exponential of a matrix, 

twenty-five years later. SIAM Review pp. 3-30. 
Rebonato, R. and M. Joshi (2001). A joint empirical and theoretical investigation of the modes of 

deformation of swaption matrices: implications for model choice. QUARC Working paper. 

Independent consultant, level3finance.com 
E-mail address: claudioalbaneseOlevelSf inance .net 

Mathematics, Imperial Collecje London. 
E-mail address: a.osseiran@iniperial.ac.uk 




